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Abstract 

We present the first complete order a s corrections to the Wilson coefficients (at the 
| high scale) of the various versions of magnetic and chromomagnetic operators which are 

ff) ■ induced by a squark-gluino exchange. For this matching calculation, we work out the on- 

C*~) [ shell amplitudes b — > sj and b — >■ sg, both in the full and in the effective theory up to order 

a 2 s . The most difficult part of the calculation is the evaluation of the two- loop diagrams 
\ in the full theory; these can be split into two classes: a) diagrams with one gluino and a 

virtual gluon; b) diagrams with two gluinos or with one gluino and a four-squark vertex. 
Accordingly, the Wilson coefficients can be split into a part a) and a part b). While 
part b) of the Wilson coefficients is presented in this paper for the first time, part a) was 
given in 1]. We checked their results for the coefficients of the magnetic operators and 
found perfect agreement. Moreover, we work out the renormalization procedure in great 
detail. 

Our results for the complete next-to-leading order Wilson coefficients are fully analytic, 
but far too long to be printed. We therefore publish them in the form of a C++ pro- 
gram. They constitute a crucial building block for the phenomenological next-to-leading 
logarithmic analysis of the branching ratio B — ¥ Xgj in a supersymmetric model beyond 
minimal flavor violation. 
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1 Introduction 



Among the rare B meson decays those induced by radiative and electroweak penguin dia- 
grams are of particular interest. Such flavor changing neutral current (FCNC) processes offer 
high sensitivity to new physics (NP) through potential new degrees of freedom beyond the 
Standard Model (SM). Additional contributions to the decay rate, in which SM particles in 
the loops are replaced by new particles such as the supersymmetric charginos or gluinos are 
not suppressed by the loop factor a/Air relative to the SM contribution. Thus, FCNC decays 
provide information about the SM and its extensions via virtual effects to scales presently 
not accessible otherwise. This approach is complementary to the direct production of new 
particles at collider experiments (for reviews see [2JE]). 

At the end of the first generation of the B factories at KEK (Belle experiment at the 
KEKB e+e" collider) [3] and at SLAC (BaBar experiment at the PEP-II e + e~ collider) [5], all 
present measurements in flavor physics including those from the Tevatron B physics programs 
(CDF [6J and DO [7] experiments) have not observed any unambiguous sign of new physics, 
in particular no 0(1) NP effects in any FCNC process [8,9|. This implies the famous flavor 
problem, namely why FCNC processes are suppressed. It has to be solved in any viable new 
physics model. It is well-known that the hypothesis of minimal flavor violation (MFV) |10pil| . 
i.e. that the NP model has no flavor structures beyond the Yukawa couplings, solves the 
problem formally. However, new flavor structures beyond the Yukawa couplings are still 
compatible with the present data [12] because the flavor sector has been tested only at the 
10% level, especially in the b — > s transitions. 

Among the penguin modes, the inclusive decay B — > X s j is the most important one, be- 
cause it is theoretically well understood and at the same time it has been measured extensively 
at the B factories. While non-perturbative corrections to this decay mode are subleading and 
were recently estimated to be well below 10% [33], perturbative QCD corrections are the most 
important corrections. Within a global effort, a perturbative QCD calculation to the next-to- 
next-to-leading-logarithmic (NNLL) level has been performed and has led to the first NNLL 
prediction of the B — > X s ^f branching fraction [14,15j which also includes the nonperturbative 
contributions. Using the photon energy cut Eq = 1.6 GeV, the branching ratio reads 

B(B -)• A s7 ) NNLL = (3.15 ± 0.23) x 10" 4 . (1) 

The combined experimental data according to the Heavy Flavor Averaging Group (HFAG) [16 J 
leads to 

B{B X sl ) = (3.55 ± 0.24 ± 0.09) x 10 -4 , (2) 

where the first error is combined statistical and systematic, and the second is due to the 
extrapolation in the photon energy. Thus, the SM prediction and the experimental average 
are consistent at the 1.2cx level. This implies very stringent constraints on NP models like the 
bound on the charged Higgs mass in the two Higgs-doublet model (II) [171118] (M H + > 295GeV 
at 95% CL) [14J or the bound on the inverse compactification radius of the minimal universal 
extra dimension model (mACD) (1/R> 600GeV at 95% CL) [19]. In both cases the bounds 



1 



are much stronger than the ones derived from other measurements. Constraints within various 
supersymmetric extensions are analyzed in [2"Urf3"T] (for overviews see Bounds on the 

Little Higgs Model with T-parity have also been presented [36]. Finally, model- independent 
analyses in the effective field theory approach without [37] and with the assumption of minimal 
flavor violation [lO^llj also show the strong constraining power of the B — > X s ^f branching 
fraction. 

While within the perturbative SM calculation NNLL precision has been achieved [14, 
1 15 j . also within supersymmetric theories higher order calculations have been pushed forward 
in recent years. In particular, the complete NLL QCD calculation of the MSSM with the 
additional assumption of minimal flavor violation was presented in [28] including a published 
computer code [38J. Beyond minimal flavor violation, the most important role is played by 
the non-diagonal gluino-quark-squark vertex due to the large strong coupling which comes 
with this vertex. This flavor non-diagonal vertex is induced by squark-mixing to the extent 
as it is misaligned with quark mixing. It represents a new flavor structure beyond the SM 
Yukawa couplings. A complete LL analysis of the corresponding gluino contribution to the 
inclusive decay rate was presented in [23]. Some "beyond LL effects" were estimated [27], but 
a general NLL analysis of the gluino contribution is still missing. 

Besides these NLL contributions due to the gluino vertex, there are of course more NLL 
corrections with non-minimal flavor violation; they involve electroweak (gaugino and higgsino) 
vertices. However, such contributions are in general suppressed compared to the corrections 
computed here. There are two types of such contributions at the NLL level: First, there are 
electroweak corrections to the non-minimal LL gluino contribution (in which the electroweak 
vertex is flavor-diagonal or MFV-like) which are naturally suppressed due to the smaller 
coupling constants and due to the CKM hierarchy. Second, there is also non-minimal flavor 
violation via squark-mixing in the electroweak vertices possible. But such contributions are 
already suppressed at the LL level compared to the gluino contribution due to the smaller 
coupling constant, apart from the chargino contributions in specific parts of the parameter 
space in which for example the trilinear coupling A^ ?J is very large. These features do not 
change of course when gluon- and gluino-induced NLL corrections are added to such LL 
contributions Summing up, the complete NLL corrections induced by the gluino vertex 
represent the dominant contribution beyond MFV at this order in most parts of the MSSM 
parameter space. They are complementary to the MFV contributions at the NLL level which 
are given in |28j. 

In the present paper we work out for the first time the complete corrections of order a s 
to the Wilson coefficients (at the matching scale fiyy) of the various versions of magnetic 
and chromomagnetic operators which are induced by a squark-gluino loop. In our procedure 
all the appearing heavy particles (which are the gluino, the squarks and the top quark) 
are simultaneously integrated out at the high scale, which we call fj>w- For this matching 
calculation, we work out the on-shell amplitudes b — > S7 and b — > sg, both in the full and 
in the effective theory up to order a 2 s . The most difficult part of the calculation is the 
evaluation of the two-loop diagrams in the full theory; these can be split into two classes: 
a) diagrams with one gluino and a virtual gluon; b) diagrams with two gluinos or with one 
gluino and a squark-loop. Accordingly, the Wilson coefficients can be split into a part a) and 
a part b). While part b) of the Wilson coefficients is given in the present work for the first 

1 Still, the leading chirally enhanced corrections can be easily calculated by inserting the effective Feynman 
rules of [32] into the results of [24] , 
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time, part a) was published by Bobeth et al. pQ. We explicitly checked their results for the 
Wilson coefficients of the magnetic operators, which were obtained using an off-shell matching 
procedure, and found perfect agreement. 

Moreover, we work out all steps of the NLO renormalization procedure in great detail; to 
our knowledge, the complete set of renormalization constants, which we give in Appendices iBl 
to El have not been presented before. We also discuss how the recently discussed chirally- 
enhanced terms |31] enter our scheme. Our results for the complete NLO Wilson coefficients 
are fully analytic, but far too long to be printed. We therefore publish them in the form of a 
C++ program. 

The remainder of the paper is organized as follows: In the next section we present de- 
tails on the framework used for our calculation. In particular, we list all relevant operators 
contributing to B — > X s ^. In Section [31 we describe in detail the calculation of the Wilson 
coefficients of the magnetic operators which is extended to the chromomagnetic case in Sec- 
tion^ All our results are implemented in a computer code which is described in Section [5] and 
we present our conclusions in Section In Appendix we present explicit analytic expres- 
sions for the LO result and Appendices [B] to [F] provide details on the NLO renormalization 
procedure. 

2 Preliminaries 

2.1 Squark-quark-gluino interaction 

The minimal supersymmetric SM (MSSM) allows for generic new sources of flavor violation 
beyond the Yukawa structure in the SM, thus, beyond the minimal flavor violation hypothe- 
sis. Besides the usual quark mixing also the squarks induce flavor mixing due to a possible 
misalignment of quarks and squarks in flavor space. 

More explicitly, in the MSSM without soft SUSY-breaking, flavor-violation is induced by 
the superpotential terms containing the Yukawa couplings of the matter superfields to the 
Higgs superfields 

£ Yukawa = -.f%(h d )ij Q°U d + U?{K)ij Q^u ~ »H d H u . (3) 

The index "o" indicates that the current eigenstate basis is used. Following the notation 
of [39U40] . we define u £ as the scalar component and as the fermionic component of the 
superfield U° (and for D° analogously); the scalar and fermionic components of the superfield 
Q° are defined as Q° = {u° L d° l ) and Q = {u° L d° L ), respectively. H u ,d are the Higgs superfields. 
The matrices h d and h u act on flavor space. Diagonalizing these matrices via bi-unitary 
transformations defines the super-CKM basis. 

In addition, there are soft SUSY-breaking terms leading to flavor transitions of the quarks. 
When restricted to the terms relevant for squark masses and quark-flavor transitions, the soft 
part of the Lagrangian can be expressed in terms of the scalar fields in the current eigenstate 
basis: 

£soft _ _ Q° m ~ _Q° — u° m rr? U i jU°^ — dftim^ .Aftj 

- (^AfjHdQidnj + A^H u Q°u% + h.c.) , 4 

where u°, d° and Q° are the squark fields and m„, and niQ the corresponding hermitian 
mass matrices. H u and H d are Higgs doublets and A u and A d are 3x3 matrices in flavor 
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space. For the latter no proportionality to the Yukawa couplings is assumed; they are left 
completely general. 

After the diagonalization of the Yukawa couplings in ([3j), i.e. in the super-CKM basis, the 
squark mass matrices are of the form (q = u or q = d) 



K = I Zl' LL Zl' LR ) + ( F * LL + D *> LL „ \ LR n ) , (5) 



Q V m q,RL m q,RR J V F 1,RL F q,RR + D q,RR 

where the F-terms from the superpotential and the D-terms from the gauge sector are di- 
agonal 3x3 submatrices of the 6x6 mass matrices M q . This is in general not true for 
the additional terms, originating from the soft part of the Lagrangian, given in Q (for more 
details see Chapter III of [H]). Thus, the (hermitian) squark mass terms 

(4 4)^ (J) and (44)^l(^)> (6) 

where the fields ul,r arid di,R now refer to the super-CKM basis, are not diagonal in general. 
Only when performing an additional appropriate unitary field transformation defined through 

£)-(&)•■ 

with 6x3 matrices T(u,d)(l,r)^ one obtains a diagonal matrix 

(T QL ,T QR )M 2 q [ J" ] , (8) 
and the squark mass terms have the form 

q\r QL ,T QR )M 2 q [~J L )q. (9) 





q represents the six squark fields in the mass-eigenstate basis. Note that in (|8|) and ([9]) one 
either has (q,Q) = (u,U) or (q,Q) = (d,D). The V matrices satisfy the following unitarity 
relations: 

E(r&rg2 + rg*r&) = s kh , fl t qx t qy = s ij s XY , (10) 

8=1 k=l 

where Q = U, D and X, Y = L, R. 

The squark-quark-gluino coupling is flavor-diagonal in the super-CKM basis; thus, the 
last rotation of the squark fields of ([9]) gives rise to the following squark-quark-gluino vertices, 
which are non-diagonal in flavor-space 

9,a mmm <^ - ig s V2Tp a (Tq L P l - Tq R P r ^ , (11) 

a, a 

9, a ■mmm J <^' - ig s V%T£p [Tq* l P r - Tq R P l ^ , (12) 
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Further Feynman rules may be found in |41H43j . 

At this point, a comment on the super-CKM basis is appropriate. While in [27] and [28J 
a definition of the super-CKM basis is used in which the physical quark mass matrices are 
diagonal, we will apply a minimal renormalization scheme to the squark mixing matrices V 
appearing in (jlip and (|12p . This corresponds to a definition of the super-CKM basis in which 
the Yukawa couplings remain diagonal also at the loop level [3TJE3JE1] • 

In the following calculation, we will express our results for the Wilson coefficients in terms 
of the diagonal squark masses rriq k and the V matrices, renormalized in the DR scheme. 

2.2 Operator basis 

In supersymmetric models, various combinations of the gluino-quark-squark vertex lead to 
|A(B)| = | A (5)1 = 1 magnetic and chromomagnetic operators (of 07-type, Os-type) with an 
explicit factor a s , and to four-quark operators, with a factor a 2 when integrating out squarks 
and gluinos [23]. These squark/gluino induced operators define the effective Hamiltonian H^ s , 
which is given in (|13j) below. 

Using the standard effective field theory approach, the aim is to resum the following terms: 

LO: a s {a s L) N 

NLO: a s a s (a s L) N , (N = 0, 1, ...), 

respectively at the leading and next-to-leading order where L denotes large logarithms of the 
ratio m 2 b /M 2 usy . 

As discussed in [55] . %e// * s unambiguous, but it is a matter of convention whether the 
a s factors, peculiar to the gluino exchange, should be put into the definition of operators or 
into the Wilson coefficients. It is convenient to distribute the factors of a s between operators 
and Wilson coefficients in such a way that the matching calculation and the evolution down 
to the low scale of the Wilson coefficients are organized exactly in the same way as in the 
SM such that the anomalous-dimension matrix indeed has the canonical expansion in a s and 
starts with a term proportional to a\. 

The effective Hamiltonian T-L^ s is of the form 

?4 = £ CisMOisQt) + c Uti°rM ( 13 ) 

i i q 

In the second term, we sum over all light quark flavors q = u,d,c,s,b. 
At dimension five the following two-quark operators contribute 

1M = eg 2 s (fi) (sa^P R b) , O' = eg 2 (v) (sa^P L b) , 

(14) 

8§rg = g s {ji) g 2 (n) {sa^T a P R b) G% , 0' 8 - g - g = g s {^) g 2 s (») (sa^T a P L b) G% , 

with Pl = |(1 — 75) and Pr = |(1 + 75). In the SUSY contributions corresponding to these 
operators, the chirality-violating parameter is the gluino mass rrig which is included in the 
corresponding Wilson coefficients: 
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At dimension six the two-quark operators with chirality violation coming from the b- or 
c-quark mass read: 



7 b,g = eg 2 s {pL)m h {n) (sa^P R b) , 0' 7h ~ g = eg 2 s (^m b (fi) (sa^P L b) , 

%Kg =g s ^)gl^)m b ^)(- S a^T a P R b)Gl v , 0'= g s {») g 2 ^)m b (p) (sa^T a P L b) G* 



[IV 1 



O 



7c,g 



egj(fi)m c (fi) {sa^P R b) F„„ , 



0' 7c - g = e g 2 ( M ) m c (») (sa^P L b) F^ , 



0^- g =g s ^)g 2 M^){- S <y^T a P R b)Gl v , Q' Sc - g = g,(jj) g 2 s ^)m c ^) (sa^T a P L b) G° 



[IV 



The following four-quark (axial-)vector operators contribute: 



(15) 



Ou,~ 9 = gKv)(^P L b)(q^P L q), 
OU, = gt{li){s al „P L bp){q^P L q a ) , 
gM(s lfl P L b) (q^P R q), 
gi(v)(saJtiPLbp) {qpl^Piiqa) , 



OIL = gM(si„P R b)(q^P R q). 



12,9 

o q ~ 

13,9 
^14,9 



o ql 

O q ' 

13,9 
^14,9 



9i(fJ-)(s a 7^PRbf3) {q^P R q a ) , 
g* s (fi)(s^P R b) (q^P L q), 



(16) 

where color indices are omitted for color-singlet currents. These operators arise from box dia- 
grams through the exchange of two gluinos and from penguin diagrams through the exchange 
of a gluino and a gluon. 

In addition there are also the following four-quark (pseudo-)scalar operators which are 
induced at leading order by box diagrams with two gluino exchanges 



O q 

U 15,9 
^6,9 
^7,9 

o q 

U 18,9 

o q 
o q 

U 20,g 



gt(n)(sP R b) (qP R q) , 
9 A s {ij){saPRbfi) {q~f)PRq a ) 
gt(n)(sP R b) (qP L q) , 
9i(^)(s a PRb/3) {qpPhqa) ■ 

9 A s{n)(qaPRbp) (spPRqa) 

9 M(qP R b) (sP R q) , 



U l5,g 

I6,g 
(0 q> 

O q ' 

U l9,g 
O q ' 

U 20,g 



gi(n)(sP L b) (qP L q) , 

gi(v)(saPLbf3) (qpPLq a ) , 

9 M(sP L b) (qP R q), 
giitAisaPhbp) (qpP R q a ) , 
9 A s {ij)(qaPLbp) (spP L q a ) , 
gM(qP L b) (SP L q). 



(17) 



These operators are shown to mix at one-loop into the magnetic and chromomagnetic oper- 
ators [23] and, thus, contribute to a LL analysis. Note that we defined Of g - and 0\ Q - (and 



also O q 'n 7., O q n t.) differently from [23] . Denoting the operators from [23] by 0\§~ and O209' 



'19,9' ^20,9 

the transformation between the operator bases reads: 



O q 

U 19,9 



\o q , rg -\d q 9 - and O q - 



-6 q - -6 q 

2^16,9 8^20,9' 



(18) 



Correspondingly we get for the Wilson coefficients: 



^15,9 - ^15,9 ~ 4C 19,9 ' 



r~<1 C<1 _ AC<1 

^16,9 ~ °16,9 °20,9 ' 



°19,9 



~ 8C 19,g ' 



°20,9 



~ 8C 20,g- 



(19) 



Analogous relations are valid for the primed operators. 
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Figure 1: Gluon corrections: The solid, dashed, curly, and solid-curly lines represent quarks, 
squarks, gluons, and gluinos, respectively. The crosses denote the possible locations where 
the photon is emitted. 




Figure 2: Gluino corrections: The crosses mark the places where a photon can be emitted. 

3 Calculation of the Wilson coefficients of the magnetic 
operators 

The various Wilson coefficients Ci appearing in the effective Hamiltonian 7i^ s can be expanded 
as follows: 

C^) = cf\») + ^C? + 0(al). (20) 

In order to extract the Wilson coefficients Ci g ^ C^- g Cj c ^ g of the magnetic operators 
Ojg^g, C>7b,g @7c,g in ([HI) and (fT5|) . we calculate the on-shell decay amplitude b — > sj in 
the full theory and match the result with the corresponding one obtained by the effective 

-® TOTOffiT* ®- 

Figure 3: Diagrams containing a four-squark vertex. 
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Hamiltonian % $ in (|13p . We require the s quark, which we treat as a massless particle, to 
be left-handed, so that only the unprimed operators are involved. The Wilson coefficients of 
the primed operators (corresponding to a right-handed s quark) can be obtained at the very 
end by interchanging L and R in the results for the unprimed ones. 

To get the LO pieces of the Wilson coefficients, we have to calculate the amplitude A(b — > 
S7) up to order a s . The corresponding result in the full theory (expanded up to second order 
in the external momenta and light masses) is given in Appendix A. In the effective theory, 
only the matrix elements associated with the magnetic operators contribute at the lowest 
order of a s . Consequently, we get for the leading-order Wilson coefficients 



r»(°) _ k-(o) 

u 7g,g - ^79,9 



r (o) _ tAO) 

£ ^ ' U 7b, g ~ «7b,g 



e-S-0 



Ci% = 0, (21) 



where the functions K^~~ and are given in Appendix A. They represent the d-dimen- 

sional versions of the LO Wilson coefficients. 

We now discuss the calculation of the NLO pieces C^~~, C^ -, C^J- of these Wilson 
coefficients. There are two types of diagrams contributing to b — >• 57 in the full theory at 
0{al): 

Full a) Diagrams containing one virtual gluon and one gluino, shown in Figure [TJ We call 
these contributions "gluon corrections". 

Full b) Contributions with two gluinos or diagrams containing squark tadpoles, depicted in 
Figures [2] and [3l We refer to these contributions as "gluino corrections" . 

As we will see, these two contributions can be renormalized separately in the full theory. It is 
therefore convenient to split also the renormalized result obtained in the effective theory into 
two contributions. This leads to a decomposition of the NLO pieces of the Wilson coefficients 
of the magnetic operators into two pieces, as e.g. C^- - = C^~'~ + ~ . 

The mentioned two contributions in the effective theory can be characterized as follows: 

Effective a) Renormalized one-loop matrix elements (s7|C- ] Oi\b) associated with the Oj- 
and 08-type operators; one-loop matrix elements of the operators 0\ x - to 0\ A - (which 
are finite) weighted with the penguin part of the corresponding LO Wilson coefficients 
(see (31) in [23]); tree-level matrix elements associated with the operators 
^C^f (s 7 \O n , g \b),i = b,g. 

Effective b) One-loop matrix elements of the operators 0\ x ~ to 0\ A ~ (which are finite) 
weighted with the box part of the corresponding LO Wilson coefficients (see (32) in [23J); 
renormalized one- loop matrix elements of Cf 5 ~Of 5 ~ to C| -O^. (see (33) in [23]); 

tree-level matrix elements associated with the operators ^ C^~ b (57107^16), i = b,c, g. 

For the precise definition of these two parts, see (f3"0j) and (f5Tj) . 

In the following, we separately work out C^' a and C^ ,b using dimensional regularization 
(DREG) both on the full theory and the effective theory side. Later we will take into account 
the important point, that one should use dimensional reduction (DRED) in the full theory 
to preserve supersymmetry. This is done in Section 13.41 by appropriately shifting the strong 
coupling constant g s and the gluino mass in the leading order contributions of the full theory. 
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3.1 Calculation of C^~ a and 
3.1.1 Full theory 

To get the order a 2 corrections to the decay amplitude corresponding to the contributions of 
class a) in the full theory, we expand the diagrams shown in Figure [1] in inverse powers of the 
heavy mass M (M = m g or M = rriq), using the hard-mass procedure [3]. We systematically 
retain all terms up to order 1 /M 2 . When evaluating a genuine two- loop diagram in Figure [1] 
according to this method, we get two contributions: The so-called hard contribution which 
is obtained by a naive Taylor expansion of the two-loop integral with respect to the external 
momenta, and the so-called soft contribution (for details how to get this contribution, see 



Adding the hard contributions of all diagrams defines the hard part iA^^u h, while the sum 
of the soft contributions defines the part %A a f u ii,si so that 

Note that in this section the symbol A is used for the contribution of order a 2 of the relevant 
amplitude. Instead of calculating the soft part iA^^u s according to the rules of hard-mass 
procedure, it can be generated using the effective theory at LO (as we checked explicitly) in 
the following way: 

i<Tull,s =4Z (si\0 7b ,g\bh-lo OP + K\% { S1 \0 Tgrg \b) Uoov 

+ (sTiOsftjI&Uop + Kfl- (vy\O m \b)^ 

21 (23) 



{Z^ a -l)K^-{ S1 \0 7Kg \b) tI 



The various (ci-dimensional Wilson coefficient) functions K^°\ which in the limit e — > 
coincide with the correponding ordinary LO Wilson coefficents , are given in Appendix A. 
Cf~ denote the penguin parts of the Wilson coefficients C?- (i = 11 — 14); they are explicitly 
given in (31) of [23]. Note that it is important to keep the functions in d = 4 — 2e 
dimensions, because they are multiplied with one-loop matrix elements which are singular in e. 
On the other hand, the matrix elements of the operators 0\ - (i = 11 — 14) are finite, therefore 
it is sufficient to weight them with the corresponding (4-dimensional) Wilson coefficients Cf'~. 

Finally, Zm^' a denotes the renormalization constant of due to gluonic corrections in the 
on-shell scheme: 



7 OS,a t a s n 
mb = 47T 



31n-^ +4 

m b 



(24) 



with C F = 4/3. 

The renormalized amplitude iA r ^ ull in the full theory is obtained by adding the counter- 
term iA^j ull induced by renormalizing the following quantities in the leading order expres- 
sion for the amplitude given in Appendix A: The gluino-mass m g , the down- type squark 
masses mj fe and the strong coupling constant of Yukawa type g s ,Y- Furthermore, the on-shell 

renormalization constants y Z®^' 11 and y Z^f' 0, have to be attached. As we are discussing 
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gluon-corrections in this section, we only take into account the gluonic contributions of the 
involved renormalization constants. The counterterm amplitude then reads 



1 A a,full 



+ 4%( S1 \o 7 - grg \b) tlee ) + m ih 



(25) 



Zg sY , derived in the Appendix IF. 1[ reads 



A 9s,Y 



1-^ 

47T 



3(C A + C F ) 
2e 



(26) 



with Ca = 3, while the on-shell renormalization constants for the b quark and the massless s 
quark, which will drop out in the matching equation (p?2"j) . read (see Appendix [Bl ([55]) ) 



Z. 



OS,a 
2b 



1 



4tt 



C F 



- + 4 + 3 In ^ 



(27) 



>yOS,a ^ 

A 2s - U - 



In (j25j) . the contribution R a _~ is obtained by shifting (according to the MS scheme) 

9^k 



m 



a 

2 2 

m j — >■ m,- 



1 - — C A - 
47r e 

i-2:o,» 

47r e 



(28) 



in the leading order expression for iAf u u in Appendix A, followed by an expansion in a s and 
keeping the term proportional to a 2 s . 

Finally, putting everything together, the renormalized result in the full theory reads 



%A 



a, full 



-4 A 



+ 



bare 
a,full,h 



+ iA 



bare 
a, full,, 



25Zl Y + 1 -6Z° b s > a + 1 -5Z os - 



J 2x 



(29) 



x (k® s (s 7 \0 7b , g \b) tica + (s 7 \0 7 - g , s \b) tiee ) + EZ 
We now turn to the effective theory side. 



gdk 
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3.1.2 Effective theory 

After converting the running 6-quark mass present in the operator 7bt § to the pole mass, the 
renormalized amplitude iA^ff can De written as 

iAg ff =g C% g a ( S1 \O rgrg \b) Uee + g C$ g a ( S1 \0 7brg \b) trcc 

+ (sjpr^lb)^ + k\% ( S7 |073,3|fc)i-ioo P 



Ec5 a ( S7 |^-|6}, 100! , 



i=ll 



(30) 



+ 



2 5Z 9s + 8 Z 77f)g + - <5Z° 5 ' a + - 8Z° S ' a 
2 5Z 9s + 8 Z 77M + - 8Z° s ' a + ^ (5Z 2 °/' a 



4&<*7|07§,§|&> t 



+ 



where the terms in the lines 2-5 correspond to i A^J^ s in (|23|) and will cancel in the matching 
equation ([32]). Those renormalization constants, whose effects do not drop out in the matching 
step (|32|) . read 



a s 23 a s 9 



a s 13 a s 16 

Z 77,6 9 = 1 + ^ — , ^87,99 = Z 87 , bg = 1 " ^ ^ ' 

Note that in ([3U|) we put the (i-dimensional versions K 7g g , ^7bg> ^-ig g anc ^ -^86 <? °^ ^he ^0 
Wilson coefficients, which implies that the infrared structure on the full and effective theory 
side are identical. Moreover, the use of the d-dimensional versions of the Wilson coefficients 
in the effective theory also implies that all contributions due to the on-shell renormalization 
constants of the b and s quark will cancel out in the matching condition (|32p . These features 
manifestly reflect the fact that the determination of the short- distance Wilson coefficients is 
independent of the infrared structure and of the choice of the external states. 

3.1.3 Extracting dfff and C$£ 

The coefficients C^~'~ a and C^ g a are obtained by requiring that 

iA r a e j u u = iA r a ™ff ■ (32) 
Using Xd k = m?~ /m~ and = ln(^ 2 /m|), we get 

n {l),a _ 4 _i_ 1 r^ r tt r b..(i),«/ \ 

" 3 16^ ™ ^ 1 drL dl h ^9 ^ ] ' 

9 k=i 

r {l),a _ 4 1 1 ST v kb v ks* Al),a ( , 
U 7b,g ~ -3 m 2 DL n 7b ,g \ X d h )> 



(33) 



m 9 fc=l 
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with 



,(i)a,, f 55x 2 + 69x + 2 2x (9x 2 + 40x + 14) \ x (x 2 - lis + 13) , 2 

f x ) = + i - '- ki L + — ; u In 2 x 

7 ^ V ' \ 9(x-l) 3 9(x-l) 4 j 18(x-l) 4 

19x 2 + 60x-15 T . , 28x 3 + 120x 2 - 15s - 1 , 

-Li 2 (l - x) H — —7 mx 



+ 



9(x-l) 3 v ; 9(x-l) 

-14x 2 - 333x + 83 



18(x - 1) 



3 



,m a , . f -5x 3 + 384x 2 + 609x + 20 x (l8x 2 + 107x + 43) , \ x(x + 11) , 2 

s (X) = ; r-. : In X L,, H ; —r 111 X 

7b > 9 V ' \ 108(x - l) 4 18(x - l) 5 J 36(x - l) 4 

-17x 2 -86x + 15 T , . 87x 4 - 537x 3 - 2997x 2 + 265x + 14 , 

Li2(l — x) H — — — -= lnx 



18(x-l) 4 v ' 324(x-l) 5 

-799x 3 + 1719x 2 + 10431x - 1847 



972(x - l) 4 

(34) 

Our results are in full agreement with those of Bobeth et al. [1] who matched the corre- 
sponding off-shell Greens functions. 

Finally, we emphasize that one could also insist to use the four-dimensional version of the 
Wilson coefficients in the effective theory, because one can work out the explicit infrared (IR) 
structure of the amplitudes in the effective and full theory by distinguishing IR poles (1/eir) 
from ultraviolet (UV) poles (1/e). One then finds that the IR structures in iA r ^ ull and iA a e ™^ 
are of the form 

4 + - ) 4 0) (0 7 ) trcc and ( 4 + - ) Cf (O 7 ) troc , 

£ ir £ ir / \ e ir £lr / 

respectively. In both cases the same singular factor multiplies the lowest order amplitude of 
the full theory and the effective theory, respectively. On general grounds, the corresponding 
statement about the IR structure holds for the bremsstrahlung contributions. So the correct 
matching conditions, to extract NLO pieces of the Wilson coefficients C^~'~ and C^'-* , can 
be achieved in this case by discarding these explicitly identified IR sensitive terms in iA a e ^ ull 
andU-^. 



3.2 Calculation of Cf>~, C$/, and cfy b 
3.2.1 Full theory 

In order to get the full theory NLO contributions to the decay amplitude belonging to class 
b), we have to calculate irreducible two- loop diagrams with two virtual gluinos depicted in 
Figure [2] and contributions with squark tadpoles shown in Figure We denote the contribu- 
tions of these diagrams by iA\ a '[^ red j ull . Additionally, there are reducible diagrams shown in 
Figure H] leading to a contribution to the decay amplitude denoted iA^^ d j ull . 

First, we discuss the calculation of the irreducible diagrams. We expand the corresponding 
two-loop diagrams according to the rules of the hard-mass procedure [44J in inverse powers 
of the heavy masses, which are in our application rriq k , rrig and mt- The class of diagrams 
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•<8K. 




Figure 4: One-particle reducible diagrams, which are taken into account through the quark- 
field renormalization constants 



containing squark tadpoles (Figure[3]) does not contain any light propagators and can therefore 
just be naively Taylor-expanded in the external momenta. The diagrams with two virtual 
gluinos (Figure [2|) can also contain propagators of light quarks. Thus, besides the naively 
expanded two-loop diagrams, we also have to take into account contributions with one-loop 
sub-diagrams containing all heavy lines which have to be expanded in the external momenta. 
The result is inserted in the second loop which contains only light propagators as an effective 
vertex. Logarithms of the small masses are only generated in the latter contributions and 
cancel in the matching procedure against equal terms on the effective side stemming from 
matrix elements of the four- Fermi operators 0\ x ~ to ~ . 

The integrals which result from the naive Taylor expansion are two-loop vacuum inte- 
grals. Using integration-by-parts identities |51U52j and partial fraction decompositions, we 
can reduce all these integrals to the following master integrals: 



d d k 



I 



(2ir) d k 2 - 
d d k x d d k 2 



m- 



(4tt)2 



' r(i-^,»"- 2 



111-, 



-4e r(l + C ) 



{2ir) d {2ir) d {k 2 - m\) {k 2 - m 2 ) [(h + k 2 ) 2 - m 2 ) 512vr 4 3 (1 - e)(l - 2e) 



1 + x + y 2(xlnx + y\ny) 



+ 



x In 2 x — y In 2 y + (1 — x — y) In x In y — \ 2 <fr(x, y) 



(35) 



where d = 4 — 2e, A 2 = (1 — x — y) 2 — 4xy, x = mf/m 2 , and y = m 2 /m 2 . The function $(x, y) 
has been evaluated in [50J. For A 2 > it is represented by dilogarithms as 



$(x,y) = - 



, 1 + x — y — A , 1 — x + y — A 
2 In ln- 



2Li. 



1 + x — y — A 



2 
1 



In x In y 



2iy * + y-\*_ 

2 3 



(36) 



whereas for A 2 < it can be written in terms of Clausen's functions 

2 



r-i / •> l + x + y \ l+x-y 

( 1-2 j Jarccos — - — — — | • ( l L > I 2arccos — j-^= — 



+ CI2 ( 2 arccos 



2^/xy 
1 - x + y 

2Vv 



(37) 



We now turn to the reducible diagrams. Their contribution to the decay amplitude can be 
obtained by amputating the external fermion legs (with squark-gluino self-energy insertions) 
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followed by attaching the LSZ-factor, i.e. the gluino parts of the renormalization constants 
of the b- and s-quark fields in the on-shell scheme which are derived in detail in Appendix IB] 
This amounts to replace the T^/i? matrices in the LO amplitude, given in Appendix lAl 
according to 

v ki . r kj 1 c 7 b,L/R _ r kj 1 f,yb,L/R,H r 7 b,L/R,AH\ , q s>\ 

1 DL/R -»■ 1 DL/R 2 6Z ji ~ 1 DL/R 2 V^H + ) ' ^> 

where a summation over the index j is understood. 

In (|38|) , we split the renormalization constants into a hermitian and an antihermitian part 

as 

6z b,L/R = 6z b,L/R,H + 5z b,L/R,AH _ (39) 
Explicitly, they read in the off-diagonal case (using Xd k = m 2 ~ /m~) 

xr7 b,L,AH «s Cf V^j/^&i -pfcj* . -r,kj* r ki \f m i +m 'jr< \ , / \) 

* = 4^ m 2 -m 2 ^ Y L D Rmj + r Di r ^ mi ) I — hM ~ ±rn~ g f 2 {x dk ) I 

+ (r&rgj (m 2 + m 2 ) + 2r^r^m,m,) / 3 (x dfc )| , 

(40) 

= ^_ Cp £ ^(rgjrg^ + r&rg^o/itoJ + r^r^/ 3 (x d j) , (41) 

where m; are the masses of the down-type quarks. The corresponding results for the right- 
handed versions are obtained by interchanging the labels R and L. The functions fi(x) are 
given in Appendix [Bj In the diagonal case, we get (using L a = ln(^ 2 /m|)) 




Performing the shifts (|38p in the LO result yields iA^™ d full a ^ ^ ne NLO level. Since the 
hermitian parts of the Z factors do not contain any inverse powers of the light quark masses, 
the corresponding contribution to i-A^Ted t u u can be expanded up to linear order in m&. The 
contributions proportional to and ml will end up in the Wilson coefficients of the magnetic 
operators of dimension five and six, respectively. 

For the contributions induced by the antihermitian parts of the Z-factors this splitting 
does not work since they also contain inverse powers of the light quark masses, i.e. chirally 
enhanced terms which have been recently discussed in |31| . Instead, we proceed in such a 
way that the contributions involving the LO functions K^~~ and will end up in the 

Wilson coefficients of the magnetic operators of dimension five and six, respectively. This 
ends the discussion for constructing iA^ a ^ d * ull . Note that in principle one could renormalize 
the r matrices in a way that the antihermitian parts of the quark field renormalization 
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constants get exactly cancelled. However, the terms induced by the antihermitian parts of 
quark field renormalization constants remain in the decay amplitude in our scheme in which 
we renormalize the T matrices minimally (see counterterms induced by the renormalization 
of the r matrices below and also Appendix IF.ip , 

We now discuss the counterterm contribution iA a ^ ull which gets induced by renormalizing 
the strong coupling constant of Yukawa- type g s ,Y, the gluino mass m g , the squark masses m dk , 
and the T D L / R matrices in the LO amplitude. We get 



1 A bjull 



--2 5Z, 



9s,Y 



K [ Z (sj\0 7b , g \b) tvcc + A ( S1 \0 7 g, § \b) trcc ) + R\; 



(43) 



bZ h n „ is derived in Appendix [FJ] and reads 



9s, Y 



sz: 



9s, Y 



4tt e 



-Cf + trrif 



(44) 



with tr = 1/2 and nt = 6 denoting the number of flavors. The term is obtained by 

replacing 

nf 



1 + 2 tr ■ 
Att e 



2 v 2 , a *ri 
m-: — > m,- H Cf 



1 



d 



(45) 
(46) 

(47) 

in the leading order expression for i A f u n in Appendix |XJ followed by an expansion in a s 
and keeping the term proportional to a 2 s . The explicit sum in (|46j) runs over all down- type 
quark flavors of mass ro^, while the one over I = 1, 6 is implicitly understood. Again, 
Xd k = rn 2 ~ /rri~ and the matrix Xrj is given by 



and 



+ —mlX%x dl X%- 



DL/R 



1 DL/R + 0L D L DL/R' 



v kk' 



rfci -pfc'j* 



1 DL A DL ! 



(48) 



with i summed over i = 1,2,3. The shift in the coupling matrices, described by the matrix 
STp 1 is derived in Appendix IF. II and reads 



5T% = -2- 



1 



71" X dl - X dk 



C F 



^kd Wd* I T-ifcd Wd* 
X ncfl npl ; 



DL 1 DR 



DR L DL. 



vkn _ x^nZ 
"gA D x dn A D - 



(49) 



for k ^= I . 



For k = I, we have <5r^ = 0. Detailed derivations of the renormalization of the squark and 
gluino masses are given in Appendices O and [Dl respectively. 

Finally, we can put everything together to form the complete contributions of class b) in 
the full theory: 



• Aren _■ Abare , ■ Abare 

l/i b,full — Ln -b,irredJuU l/i b,redjull 



+ 2 5Z] 



9s,Y 



K Z WieV^U + K^ g { S1 \0 7 ~ g , g \b) t 



'(0) 



R 



(50) 



gd k v 
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3.2.2 Effective theory 

Since the contributions of class b) do not contain any gluon corrections, the only one-loop 
contributions in the effective theory are due to the box-operators in (|17p and the box-induced 
part of the operators in (|16p (given in expression (32) in |23j). The renormalized amplitude 
iA^ e ™ff can be written as 



%A 



% C% g b ( S1 \0 7M \b) trcc + g C%> g b ( S1 \0 7brg \b) tICC + g C<»* (s^\0 7c , g \b) trac 

14 20 
i=U q i=15 q 

+ (c^'-SZxw + C$pZ 16 , n ) (sj\O rb ,- g \b) tiee (51) 



+ (C$pZ 19 ,n + C^pZ^nj {s~(\0 7Kg \b) %tee 
+ (c$pZ 19<7c + C^5Z 20Jc ) (s 7 \O rcrg \b) txee , 



with 



x „ a, 11 a, 11 

47r e 6 47r e 2 



An e 2 47r e 6 

^19,70=^-, 5Z 20 , 7c ----. 

3.2.3 Extracting cgl 6 , C^f and C«l 6 

Prom the requirement that (|50p and f)51 1) produce the same results, we can extract the next- 
to-leading order Wilson coefficients Cj~~, - and C^~. Unfortunately, the results are 
much too long to be printed. We therefore provide a computer code (see Section [5|), which 
enables the user to evaluate our results numerically for arbitrary input parameters. 

3.3 Results for Cg § , CgJ g , and 

We are now ready to recombine the contributions to the Wilson coefficients of both the gluon 
corrections of class a) and the gluino corrections of class b) as 

r W _ r (l),a r (l),b 

These results were, however, obtained by using dimensional regularization (DREG) through- 
out the calculation followed by renormalizing m g , m 2 ~ , g s y and Tdl/r in the MS scheme. 
But, since DREG introduces a mismatch between the bosonic and fermionic degrees of free- 
dom, the present calculation in the full theory has to be worked out in dimensional reduction 
(DRED) [45,46], which is expected to preserve supersymmetry. In principle this requires ei- 
ther the calculation of additional graphs featuring e-scalars or working out the Dirac algebra 
in d = 4 dimensions. 
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At leading order in a s , the results in DREG and DRED are identical, while at next-to- 
leading order the difference between DREG and DRED can be obtained by shifting the strong 
coupling g s and the gluino mass m g in the leading order contribution. Exploiting this, we 
chose in a first step to evaluate all diagrams in the full theory strictly in DREG, followed by 
renormalizing , rrig and g s in the MS scheme. In Section 13.4^ we will discuss the shifts in 
g s and m g . 

Additionally, the calculation in the full theory was performed without taking into account 
the decoupling of the heavy particles in the running of the strong coupling. This will be done 
in Section [331 



3.4 Transition to the DR scheme 



Since supersymmetry is broken in DREG, we have to distinguish between the Yukawa-type 
coupling g s y appearing in the squark-quark-gluino vertex and the gauge coupling g S} G- Only 
supersymmetry guarantees that both versions of the strong coupling are equal, which is the 
case in DRED accompanied by minimal subtraction (which is usually called DR scheme). 

The result for the NLO full theory contributions obtained in the MS scheme is expressed in 
terms of gfy an d m}~ s . The transition to the expression depending on the corresponding DR 
quantities amounts to performing the following replacements in the leading order contributions 

USES] | 



MS 

9 S ,Y 



a™ (C A - C F ) 
Air 2 



TO; 



a 



(54) 



This leads to finite shifts of the NLO Wilson coefficients which we denote by 5C^~ g ° R and 



<5Cl^~ DR . Since we treat the external b and s quarks in the on-shell scheme, their masses do 
not need to be shifted. 



3.5 Decoupling of heavy particles 



After performing the shifts in (|54p . the result on the full theory side depends on o^^- 
However, on the effective side, there are only five active quark flavors and the coupling is 



renormalized in the MS scheme. Hence, we have to express a <f ^ in terms of a^-^s on the ^ un 



„(5) 



theory side. The relation between a 



na and a fL 



(full) 



S,DR 

can be found in 



as 



a 



(full) 



(5) 
U S,DR' 



(c s 9 rr 



(55) 



with the decoupling coefficient 

.,(5) 



(err 



a 



1 



,DR 



7T 



- — In -— — TT ( In — h In ^» 

6 mf 24 \ ^ mf ^ m 2 - 

* \A:=1 u k k=l d k , 



2 Note that g s ,G does not appear in our LO expressions. 



— In ~— — TT 

2 TO? 



(56) 
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We now express or^L* in terms of according to 



(5) 



(5) (5) 

or = q v 

S,DR S,MS 



1 + 



o 



(5) 



C A 



4vr 3 



.,(5) 



Combining ([55]) and ([57]) we can finally give the relation between ct s and o^ s as 



(full) (5) 

a - = a - 

S.DB S.MS 



1 + 



a 



(5) 



7T 



^ + il n ^ + _L[yi n ^_ + yi n 

19 fi m 9/1 W-J m? ^ 



12 6 



2 1 



fe=i 



A* 

2 



2 m| 



(57) 



(58) 

This replacement leads to a finite shift of the NLO Wilson coefficients which we denote by 
6C$f c - and 5C$f c : 



3.6 



Final result for C^- g , C^-, and C) <g 



y(i) 



We now have all ingredients to give the final result for all NLO Wilson coefficients of the 
magnetic operators. They read 



(1),DR 

7b,g 

(1),DR 

7c,9 



_ n (l),a n {l),b Ar ,(l),DR Ar ,(l),dec. 



C, 

a 



_ r {l),a r (l),b Ar ,(l),DR 

_^(l),a r (l),6 
-°7c,9 + °7c,g ■ 



SC. 



{I), dec. 
7b,g 



(59) 



4 Calculation of the Wilson coefficients of the 
chromomagnetic operators 

In order to extract the NLO Wilson coefficients of the chromomagnetic operators, we have 
to match the full theory NLO result for the amplitude of b — > sg to the one in the effective 
theory. Since the Dirac structure of the photon vertices is the same as the one of the gluon 
vertices, we can use the results for the diagrams already calculated for b — > S7 and supplement 
them with the appropriate color factor. In addition to these diagrams, we have to take into 
account all those where the gluon couples to gluino and gluon lines. 

In full analogy to the procedure for extracting the Wilson coefficients of the magnetic 
operators, we also divide the calculation into the classes a) and b) introduced in Section [3] 

4.1 Results for and 

Since the calculation for the gluon corrections has been presented in [1] and we cross-checked 
the results for b —> 57, we do not repeat the calculation and directly take the result from [1]. 
It reads 



a 



(l),a = 

™ 3 16vr 2 m 



1 1 X^r kb r ks * h (1),a iV , N 



k=l 



r (X),a _ 4 1 1 \~^ r kb r ks* h (l),a, n 
°86,g - -3 7^2 2^ L DL 1 DL n 8 b,g \ 



(60) 



9 k=l 



18 



with 



7 ( 79s 2 - 12s + 5) 7s (18s 2 + 19s - 1 



12(s - l) c ' 



6(s - l) 4 



(-359s 2 + 339s - 204) 
+ ~ 24(s^TP ^i 2 (l-s) + 



ln(s) L, 



s(l81s 2 - 236s + 3l) 



48(s - l) 4 

2419s 3 - 771s 2 + 453s + 11), , x 1667s 2 - 990s + 379 

' ln(s) + - 



ln 2 (s) 



48(s - l) 4 
(747s 2 + 640s - 43) 



24(s - l) 3 



+ 



+ 



+ 



, , v 779s 3 + 7203s 2 + 93s - 11 \ r 
288(s - I) 4 j 

s 



48(s - l) 5 

45s 3 + 1208s 2 - 901s + 570) , s (45s 2 + 358s - 49) - 
— — i-Li 2 fl - x) V .„,, — '- ln 2 (s) 



96(s-l) 4 v ; 192(s-l) 4 

(-183s 4 + 30027s 3 - 10692s 2 + 6115s + 77) 



864(s - l) 5 
5359s 3 - 241425s 2 + 143253s - 59251 
5184(s - l) 4 



ln(s) 



(61) 



4.2 Calculation of C^, Cj»> b , and 

Apart from additional diagrams, there are two new features in this calculation on the full 
theory side: Firstly, the emitted gluon couples with the gauge version of the strong coupling 
g s ,Gi which, at leading order, is different from the Yukawa version. Secondly, the gluon field 
has to be renormalized. The required renormalization constants are derived in Appendices IF. 21 
and [HI respectively. They read 



sz: 



9s, G 



5ZS 



a s 
4tt 



1 C A 



a s 

- 2 ir 

4"7T 



1 C A 
- 3 tm f + — 



(62) 



Multiplying the leading order contribution on the full theory side, they both induce counter- 
terms at the next-to-leading order. However, they appear in the combination 



sz: 



+ lsz" 



0, 



(63) 



and therefore their effect cancels. 

The calculation of the decay amplitude for b — > sg on the effective side is very similar 
to the one for the b — > 57 transition. The non-vanishing renormalization constants for the 
operators read 



SZ 



15,86 



Os 11 

4vr e 2' 



SZ- 



20.86 



4vre2' 



SZ- 



20,8c 



Os 11 

4tt e 2' 



(64) 



After the extraction of the Wilson coefficients, the results are expressed in MS entities 
and we need to replace them by their DR counterparts. In contrast to the 6—7-57 case, the 
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Parameter 


C++ variable 




, LTT Ft ) 


std: : complex<double> GU[6] [6] 


(rnr, 

V UJ-j ) 




std: : complex<double> GD[6] [6] 


fm 2 , 


. . . , m 2 ~ ) 


double mupsquark[6] 


H' 


•••' m i 6 ) 


double mdownsquark [6] 


m 5 




double mg 


m t 




double mt 


m b 




double mb 


m s 




double ms 


m d 




double md 



Table 1: Input parameters. 



leading order contribution now also contains the gauge coupling g™g which has to be replaced 
according to [48j 



MS 

9s,G 



9s 



1 



o^Ca 
4tt 6 



(65) 



The decoupling of the heavy particles and the shift back to the MS scheme works exactly 
the same way as for b — )■ sj. 



5 Results for the Wilson coefficients at the matching scale \±w 

As mentioned above, the analytic results for the NLO Wilson coefficients are too long to be 
printed. We therefore provide a C++ code available at arXiv . org which enable the user to 
evaluate our results numerically for arbitrary input parameters. We also include the results 
for the LO Wilson coefficients that were given in |23j . 

The file bsg2gluino . h contains the definitions of the functions for all LO and NLO Wilson 
coefficients and can be included into any C++ codeH All input parameters, namely the 
rotation matrices of the squark fields ^m,D)(L,R) defined in Q, the squark masses, the gluino 
mass nig, the top mass mt, the bottom mass nib, the strange mass m s , and the down mass 
md, are declared as global variables within bsg2gluino .h according to the naming scheme 
shown in Table [TJ 

The names for the provided functions for the Wilson coefficients are listed in Tables [2] 
and El We define the expansion in a s as usual by 

Q = cf + g C «. (66) 

The NLO Wilson coefficients explicitly depend on the matching scale jiw Hence, these 
functions have to be provided with a value for fiyy as an argument. 



s An equivalent FORTRAN version is available from C.G. 
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Wilson coefficient 


C++ function 




r (o) 

^*77,7, 
199 


std: 


: complex<double> 


C7ggLD() 


Mo) 
199 


std: 


: complex<double> 


C7ggLDPrime() 


r (o) 

log 


std: 


: complex<double> 


C7bgLD() 


MO) 

log 


std: 


: complex<double> 


C7bgLDPrime() 


Sgg 


std: 


: complex<double> 


C8ggLD() 


MO) 

U 8gg 


std: 


: complex<double> 


C8ggLDPrime() 


MO) 
^8bg 


std: 


: complex<double> 


C8bgLD() 


MO) 
°8fcg 


std: 


: complex<double> 


C8bgLDPrime() 




std: 


: complex<double> 


C'i'g'q'LOO, e.g. CllguLOO 




std: 


: complex<double> 


C'i'g'q'LOPrimeO 



Table 2: C++ functions for the LO Wilson coefficients with i = 11, . . . , 20 and q = u,d, c, s, b. 



Wilson coefficient C++ function 





std: 


: complex<double> C 


i 'ggDRbar (double) 




std: 


: complex<double> C 


i ' ggDRbarPrime (double) 




std: 


: complex<double> C 


i'bgDRbar (double) 




std: 


: complex<double> C 


i ' bgDRbarPrime (double) 




std: 


: complex<double> C 


i ' cgDRbar (double) 




std: 


: complex<double> C 


i ' cgDRbarPrime (double) 



Table 3: C++ functions for the NLO Wilson coefficients with i = 7,8. 

To illustrate the use of bsg2gluino . h, we also provide the compilable file template . cc in 
which the header bsg2gluino . h is included and all input parameters are set to some specific 
values. The compiled program will then generate output for all included Wilson coefficients 
as shown in the file output.dat. 

6 Conclusion 

The main purpose of this paper is the two-loop matching calculation of the gluino-induced 
contributions to the processes 6—7-57 and b — > sg. While the contribution involving one 
gluino has already been available in the literature, the two-gluino part is new. We discuss 
in detail our renormalization procedure and the issues related to dimensional reduction and 
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the decoupling of heavy particles from the running of a s . Our results are presented in the 
C++ file bsg2gluino . h which comes together with this paper and allows for a convenient 
numerical evaluation. 

The results obtained in this paper constitute a crucial building block for the phenomeno- 
logical NLL analysis of the branching ratio B — > Xgj in a super symmetric model beyond 
MFV because the NLL gluino contributions are the most important corrections beyond MFV 
at this order in most parts of the parameter space in such models. 
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Appendix 



In Appendix [A] we give the LO expressions for the Wilson coefficients (in d dimensions) 
which are used at various places in this paper. Appendix [B] is devoted to the evaluation of 
the renormalization constants of the quark fields and quark masses in the on-shell scheme. 
Appendices O to [E] deal with the corresponding renormalization constants related to the 
squarks, the gluino and the gluon. For completeness we give all these quantities in the on- 
shell scheme, although for some of them only the singular piece (in e) is needed in the main 
text. Finally, in Appendix[F]we derive the 1/e pieces of the renormalization constants for the 
strong coupling constants g s y and g s ,Gi as well as for the T DL i R matrices. All the results in 
the Appendix are obtained by using DREG. 



A Leading order result for b — > S7 in the full theory 

Using the abbreviations x dk = m~ /rn~ and = ln(^ 2 /m~), the decay amplitude A(b — > 57) 
in the full theory to leading order in QCD (corresponding to the diagrams with one-loop 
gluino/squark exchange) can be written for a massless s quark as 

iA full (b -+ sj) =K$- ( S1 \0 7brg \b) tICC + k\°1- (s 7 \0 7g , g \b) tICC 

+ kZ ( S1 \0' 7b Jb) tiee + KW (s 7 \0' 7grg \b) tica . 

The first two terms on the r.h.s. of (|67p correspond to a left handed s quark and the last 
two terms to a right handed s quark. The various coefficients read (retaining terms of 
order e 1 ) 

K t = t °fif° lQd A h (-< + 6x l - ^ - 2 ) - 1 »m 

16vr 2 m| (x dk - 1) [12 2 



+e 



{Xdk 72 1} {^l ~ 49x dfc -22 + 6 (x 2 dk - 5x dk - 2) L„, 



+ Y^ x d k H x d k ) [x 2 dk - 6x dk - 6 - 6L M ) + -x dk ln 2 (x dk 



K. 



(o) 

7 9,9 



E 

k=l 



Cy ^ 1 PR 1 DlSd 

16TT 2 mg (x dk - iy 



2 i x X 



1) - x dk \n{x dk ) + e 



4 ( x d. k 



1) (3 + 2L A 



-x dk ln(x dk ) (x dk + 2 + 2L M ) + -x dk ln 2 (x dk ) 



(69) 



(0)' (0)' 

The corresponding primed coefficients Kj b - and can be obtained by interchanging 

L f-> R in the unprimed versions. Note that we explicitly kept terms of order e . For the 
individual equations ([23]) , ([29]) and ([30jl even the e 2 contributions of these coefficients would 
be required. In the matching equation (|32p . they drop out, however. 

The decay amplitude A(b — > sg) in the full theory to leading order in QCD (corresponding 
to the diagrams with one- loop gluino/squark exchange) can be written in complete analogy 
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to (|67p . The corresponding coefficients read 

6 



K. 



(o) 

8b,g 



E 

k=l 



1 DL DL 



167T 2 m? (x dfc - l) 4 



Cf - ypA 



{~ x d k + 62 4 - 3a; d fc - 2 ) - ^<4 ln 0d* 



+e (iv d * kfcdj (~ 6L ^ + x i ~ 6xd * - 6 ) + 4^4 m2 (^dJ 



72 



(x dk - 1) (6 - 5x dk - 2) L M + - 49x dfc - 22} 



-x\ \n{x dk ) + — (2x 3 dk + - 6x dk + l) 

(1 1 
-J^ x d k Hxd k ) (6L^ + 2z dfc + 9) + -x 2 dk \vL 2 (x dk ) 

+ ^ (x dfc - 1) (6 (2x 2 dk + 5x dk - 1) L M + 22:4 + 49 ^ " 5 )) } 



(70) 



(0) 



E 



1 LTR 1 DL 



\ 16n 2 mg (x dk - l) c 



Cf — 



o ( x \ - 1) ~ x d k Hxd k ) 



+e 



' ' 2 l)(2L M + 3)-^ dfc ln(x dfe )(2L At + x dfc +2) + ^x dfc ln 2 (x dfc 



+ \c A 



1 



1 



: ln(x dfc ) + - (3x dfc - 4x dfc + l) 



+ e ( ^ x d k l n 2 {x dk 



-x\ ln^J (2L M + 3) 



(71) 



+ 4 " 1) ((6a*» - 2) L M + 7x dfc - 1) 

The corresponding primed coefficients K^ b and K^- ~ can again be obtained by interchanging 
L f-> in the unprimed versions. 



B Quark-field and quark-mass renormalization in the on-shell 
scheme 

We start with the free (bare) quark Lagrangian for the down-type quarks 

C = i $ >L 7 ^<%V°z + i $U -fdArfU ~ $W<l ~ ^Idmfr^n , (72) 

and express the bare quantities (marked with the superscript (0)) in terms of renormalized 
ones: 



(73) 



m? = m,- — Srrij . 
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(a) (b) 
Figure 5: Contributions to the bare quark self-energies S^j. 

For the renormalized quark self-energies S^| n we get 

m n =S i4 - 1 -{5Z L + 6Z L %4P L - \{5Z R + 8Z R %rfP R 

+ I (5Z R % mi P L + \ (5Z L % mi P R + \ mj SZh P L + \ m, 6ZS P R ^ 
Z Z Z J Z J 

- brrii Sji Pl - 5m* Sji Pr , 

where Ejj denotes the bare quark self-energies which can be decomposed as 

E i4 = j Pl ^{P 2 ) +i/Pr ^(p 2 ) + Pl A%(p 2 ) + P R A^p 2 ) . (75) 

We note that for j ^ i the bare self-energies at order a] receive only contributions from 
diagram (b) in Figure [5] with a squark and a gluino in the loop, while in the diagonal case 
j = % also diagram (a) with an internal quark and a gluon contributes. 

B.l Off-diagonal case 

For j ^ i there are only contributions from the gluino-squark loops shown in Figure [5] (b). 
Because of this, we supplement the renormalization constants with the label b. The on-shell 
renormalization conditions read 



ReEJ? n u(pi) = 
Reu(^)E^ n = 



P?=m? ' (76) 

9 9 ' 

p A =m z 



where the symbol Re stands for taking the dispersive part. From these conditions, the wave 
function renormalization constants can be fixed uniquely. The generic formulas for the her- 
mitian and antihermitian parts of 5Z^ read (written in terms of E^, Ej|, Aj i} and A R ) 



5Z f ,H = m 2 1 _ m 2 Re - A R (m 2 )] + m^m 2 ) - A%{m 2 )] 

1 3 

+m i m j [J:f i (m 2 ) - ^(m 2 )} + mfE^mf) - m|Ej(m|)} , 
6Z f' AH =^f^ {mlA^m 2 ) + A R {m 2 )) + m^im 2 ) + A%(m])] 

1 J 

+m i m j [Y,f l (m 2 ) + Ef^rnf)] + m?E£(m?) + mfs£(mf)} . 



(77) 



The right handed versions can be obtained by replacing all labels R f-> L. Due to the unitarity 
of the r matrices, the off-diagonal renormalization constants turn out to be finite for d = 4. 
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They are needed to order e° and up to linear order in the quark masses rrii and rrij . Explicitly, 
we find 

6 r 



sz 



b,L,AH a s 



C F 



j> 



Air to? 



to' 



E a- 



DL L DR m J + 1 DL DR m i 



kj*-pki 



3 k=l 
ki -pfcj* I ^2 



to? + rrij 
rria 



fi(x dk ) - 4mgf 2 (x dk ) 



fc=i 



(78) 
(79) 



where the functions fi(x) read 

h(x) 

f2(x) 



2x ln(x) 



1)3 
.-1 

2x(x - 2) ln(x) 



x ln(x) 



(80) 



4x 



2(x - 1) 



B.2 Diagonal case 

For j = i the on-shell renormalization conditions read 

Re^[ n u{ Pl ) = 

Reu( P i)m n = 
V + mj 



9 9 



Re 



p2 _ m 2 



(81) 



Reu(pi) £■ 



+ nrii 



0. 



TO, 



We note that these conditions do not uniquely fix the renormalization constants. We are 
therefore free to choose the antihermitian parts SZ^- ,AH and SZ^' AH to vanish: 



SZ, 



L,AH 



0, 



§Z RAH = q _ 



(82) 



This is the only choice where the Srrii do not depend on the first derivative of E^, £ 



R 



A^, A^ with respect to p 2 ; in other words, the Srrii are only related to the pole position of 
the renormalized propagators, as it should be. With this choice, the other renormalization 
constants are uniquely fixed by the renormalization conditions (|8ip . 

The results for SZ^ ,H (SZ^' H is obtained by exchanging R -R- L) and for the mass shifts 
Srrii, again written in terms of £^-, Sj^, and their derivatives with respect to p 2 

(denoted by a dot) read 

(E|(m?) + Sf(m?))+^(m?)l, 



Srrii 

szl 



Re 



L 2 

Re S,v(to/ ) + rrii Re 



mi E£(m?) + rrii Eg(m?) + + ) 



(83) 
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These quantities are again needed up to order e and up to linear order in the quark mass 

As already mentioned, the diagonal quantities £|j, E^, and get contributions both 
from squark-gluino loops and from quark-gluon loops. We therefore write 5Z U ' and Srrii as 
a sum of two contributions: 

5Z^ H = 5Z^ L ' H + 5Z^ H , Sun = 5m? + 5m b , , (84) 



where the first and second term on the corresponding r.h.s. result from the quark-gluon and the 
squark-gluino loops, respectively. Note that SZf- L ' H , 5Z^ L,H , 5mf and 5m\ contain 1/e poles. 
Explicitly, we get (using the Feynman gauge for the gluon propagator and = ln(/U 2 / m|)) 



SZg** =5Z^ H = -EjZCj? (\ + f + 4 + 3 Intf/mf 



6 , x (85) 

Att~" ^ \ m- c 



k=l 



«s ~ /3 . „ , , 2 / 2- 



5m? =- t C F mi - + 4 + 3 ln(^/m 
47r V e 

1 6 

Sml =-|c F ^ (WrMi + rg fl rg&) / 3 (^J - ±m § rg £ r& / 2 (^j) (86) 
fe=i 

In (|85|) the symbol ej r has been introduced in order to mark the infrared poles and distinguish 
them from the ultraviolet ones. 

C Squark-field and squark-mass renormalization in the on- 
shell scheme 

We start with the free (bare) Lagrangian for the down-type squarks 

c = -$ t a fcfc /n& - H/4> k m, (87) 

and express the bare quantities (marked with the superscript "0" ) in terms of renormalized 
ones: 



~>S> 



K) 2 = KU 2 -K- 
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Figure 6: Contributions to the bare squark self-energies "£ k 'k- 



For the renormalized squark self-energies Eu? we get 

=t k , k - l -(5Z + 5ZA p 2 + \{8fr) k , k m\ + \m\ (5Z) k , k - 5 k < k 5m 2 -. , (89) 

Z \ / k'k Z k z fc' fc 

where T, k i k denotes the bare squark self-energies. We note that for k' ^ k the bare self-energies 
at order a\ receive contributions from diagram (b\) (see Figure [6]) with a quark and a gluino 
in the loop and from diagram (62) with a squark tadpole. In the diagonal case k' = k diagram 
(a) with an internal squark and a gluon also contributes. 



C.l Off-diagonal case 

For k' 7^ k the on-shell renormalization conditions read 



Re Ej^(p 2 = m\) = , Re £j^(p 2 = m\) J= , (90) 



where the symbol Re stands for taking the dispersive part. From these conditions and from 
we get 



SZk'k — — 2 q— Re Hfc'fc(mT ) , 

m j ~ m j 

2 ~ - (91) 

m ~ — m ~ a y 
dy d k 

When looking at the explicit result for the off-diagonal bare self-energies T, k i k below, one 
finds that these two equations are compatible with each other. Therefore only one of them is 
needed. 

The off-diagonal bare self-energies T> k i k (p 2 ) can be decomposed as 

W)=2&V) + £a ) , 02) 

where sj^(p 2 ) and correspond to diagrams (61) and (62) in Figure As indicated 

in the notation, sSl is independent of the external momentum. Explicitly, one obtains for 



28 



flj^k (p 2 ) (retaining only linear terms in the down-type quark mass m^) 



a , 



7T 



-+2 + 1,+ -^-1 In 



m~ 



2vr 



-{j? - 2m?) + 2p 2 - 3m? + (j? - 2m 2 g)L fX 



ln 1 



P 



mi 



(93) 



while reads (using Xff' = r^T^ - F^jT'fy* with i summed over i = 1, 2, 3) 



fc'fc 3vr d i 



1 1 , M 
- + 1 + In A r 

e m~ 

di 



(94) 



These expressions for and hold for both the off-diagonal and the diagonal case. The 
explicit result for SZ^'k reads (for k' / k, the label b indicates that there are only contributions 
from class b)) 



5Z b k , k = - 



x I - + 2 + L u - 



C F 



E 



m d /-rfcd*-n*;'rf 



/-pkd*-pk' d I -pkd*-pk' d \ 
\ L DR DL ' DL DR) 



m„ 



X dk ~ 1 



X d k 



In 1 1 - x dk 



1 



+^x£ , ar d| Xtf I - + l + ln^V 



1 



m 



The hermitian and antihermitian parts are: 

2 



= - - C F 



J fc'fc 



7T 



Vm"; I 1 Di? 1 DL + 1 DL 1 Di?J y 



^Jfe'Jfe 



2x dy 
2 



In 1 1 - x d , , I - 



x dk - 1 



2x, 



d k 



In 1 1 - ar d . 



7T - x dk 



X dy 



m d (T,kd*Yk'd 



. a, 
- 1 



i -pAccfe-pAi'd N 
Di? 1 DL + 1 DL 1 Di?, 



In 1 - x d - 



2x 



<4 



In 1 1 - x d 



1 



+ I - + 1 + In 



/A- 



A 4 



m- 



(95) 



(96) 



(97) 



where we used the shorthand notation x^, = m? /m? The result for (5^^ is given in the 

d k ' y 

following subsection. 
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C.2 Diagonal case 

For k' = k the on-shell renormalization conditions read 



Re£^V) = 0, for p 1 ^ mi , 



ReXr fc V)=0, for mi. 



9 2 

rr — m- 



These conditions fix 5m 2 - to be 

dk 



(98) 



(99) 



However, they only fix the hermitian part (real part) of bZkk- We choose the antihermitian 
part (imaginary part) to vanish and get 



SZ, 



kk 



dRet kk (p 2 ) 



Qrpl 



(100) 



As already mentioned, the diagonal bare squark self-energies get contributions from 
diagrams (a) and diagi 

and 5Zkk accordingly: 



diagrams (a) and diagram (b) in Figure El It turns out to be convenient to decompose 6m 2 ~ 



bZ k k = SZ kk + 5Z kk , 5m\ = 5m 2 £ + 5m z J 



2,6 
dk 



(101) 



Explicitly, we find 



W-' a =-^C F mj~ - + 7 + 3 In 



d k 47T 



2 

mi 



(102) 



dk ■ 



Smf = - ^C F 

dk 7T 



£ (m- g m d crg^r« + rg£r&) 



i+2 + L M -^-^ln|l-x d J 



a Sr , 2 



-(x dk - 2) + 2x dfc - 3 + (x dk - 2)L^ 



— m~ g X D x dl X D 



+ 1 + ln 



9 

ml 



2x dk + 1) 



In 1 1 - Xd h | 



(103) 



4 7T 



(104) 



SZkk — - 7T~Cf 

1 IT 



+ ^C F 

7T 



- + 2 + 

e 



In 1 1 - x dk | 



(2x di 



x <f fc 



1) 



. <2 



I 1 ZXR 1 DL 



x dk {f-Xdk) 
(x dk +ln|l -a d J) 
x d k 



(105) 
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' uuuuu 





U-2) 




(a) - - 

(b) 

Figure 7: Contributions to the bare gluino self-energy S. 

D Gluino-field and Gluino-mass renormalization in the on- 
shell scheme 

We start with the free (bare) Lagrangian for a gluino, which is described by a four-component 
Major ana field ip^: 

£ = \ ^II I^Il + \ $U l"d^l R - \ ^ g , R m^l L - \ ^ iL (m?) V?, r , (106) 

and express the bare quantities (marked with the superscript (0)) in terms of renormalized 
ones: 



(107) 



m~ = wig — 5m g . 

Due to the Majorana nature of the gluino, we have = (^P^lJ an d therefore 6Z R = 
{8Zi?)* . For the renormalized gluino self-energy S ren , we get 

s ren =£ - 1 (5Z^ + 5Z^'*) p 1 + m g 5Z^ P L + m g 5Z^'* P R - 5m g P L - 5m* g P R . (108) 

The bare gluino self-energy S, which gets contributions from the three diagrams in Figure 
can be decomposed as 

£ = p 1 P L £ V) + p> P R Z R (p 2 ) + P L A L (p 2 ) + P R A R (p 2 ) . (109) 

The on-shell renormalization conditions are formulated in the same way as for the quarks 
in the diagonal case j = i in Subsection IB. 21 The generic formulae for the hermitian and 
antihermitian part of Z^ and for the mass shift 5m g read 

5Z^ H = E L (m?) + 2 m? t L (mj) + m g A L (m~) + m g A R (mj) , 

^ = 2^(A*(mS)-^V|)), (110) 

5m g = m g S L (m?) + \ [A L (m 2 § ) + A R {m])] . 
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We now turn to the explicit results which we decompose as 



Srriq = <5m~ + 6m$ , 



SZi 



5Z^ L + 5Z 



7b,L 



(111) 



according to the contributions shown in Figure \7\ We find (using the Feynman gauge for the 
gluon propagator) 



a,L,H 



6Z. 



sr^a,L,AH 

5Z 3 ' 



— Ca — I- finite, 
4ir e 



5z: 



,L,H 



a 



ni 



:ii2) 



— 2 tr -± + finite, 
47T e 



dZ 3 = finite, 

where 'finite' denotes ultraviolet finite pieces not needed in this work. Numerically, we have 
tr = 1/2, Ca = 3, rif = 6. For the mass shift brxig the explicit results read in the on-shell 
scheme 



a 



<5m~ =— Ca m 



4tt 



- + 31nA,+4 

e mj 



5ml 



a s 

— 2tr — 

Air e 



nig + — tr m„ 
4ir 



QL L QL + 1 QR L QRJ 



m 



1 + ln 



m. 



m. 



+ 



Qi 



m~ 



m. 



1 + ln 



ml 



ReB (m~;mQ k ,m Qi ) 



k,i,C 



m, 



(113) 



The indices in the sums within the expression for 5m b ~ run over the following ranges: k = 
1,...,6, i = 1,2,3, and Q = D,U. For Q = D (Q = U) the symbols mQ t and denote 
the masses of the down-type (up-type) quarks and down-type (up-type) squarks, respectively. 
The function Bq appearing in (|113p is defined as 



Bo(p 2 ;m l5 m 2 ) 



1 , n — P 2 x (1 — x) + m? x + m% (1 
dx In — - ^ 



iS 



(114) 



E Gluon-field renormalization 



To fix the Wilson coefficients of the chromomagnetic operators Og, which are related to the 
process b — > sg, we need the on-shell renormalization constant (LSZ-factor) for the gluon 
field. The diagrams which define the bare gluon two-point function at order a s are shown 
in Figure O The corresponding gluon field renormalization constant Z3, defined through 
Z3 (where and A^ denote the bare and renormalized gluon field, respectively), 



— v ~s ^ 
can be written as 



Z 3 = 1 + bZl + bZ\ , 



(115) 
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(a) (b) 
Figure 8: Contributions to the bare gluino two-point function. 



according to the contributions shown in Figure [HJ As we do not want to reproduce the order 
a s correction to the Wilson coefficients Cg due to gluons (which were computed in pQ), we 
only need bZ\ in our work. We obtain 



5Z\ = 
3 4vr3 




k=l 




(116) 



where nt = 6. 

F Renormalization of g s ^y, ^dl/r an( ^ 9s,G iri the full theory 
F.l Renormalization of g sY and Tp L , R 

We use the squark-quark-gluino- vertex in order to compute the one-loop counterterms for the 
strong coupling constant g s y of Yukawa type and the T DL / R matrices in the MS scheme. At 
tree-level, this vertex for a down- type squark of flavor k (color /?), a down- type quark of flavor 
i (color a) and a gluino (color a) reads (see (fTTj) ) 

Vt\™ = -i9sV2T$ a (Y h i L P L - T% R P R ) . (117) 

At one-loop order, there are four one-particle irreducible diagrams contributing to the vertex 
correction. Three of them contain a virtual gluon in the loop, while the remaining one 
contains a quark-squark-gluino loop. It turns out that the latter contribution is finite (after 
using the unitarity relations) and thus there is no contribution in the MS scheme. For the 
vertex correction SV^ red due to these irreducible diagrams we obtain (retaining only the 
UV-singularities) 

gyjrred = ^ ^C A +C F ^ 

47r e 
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Additionally, there is a counterterm contribution 5V£f which contains the following sources: 
renormalization of the quark, squark and gluino fields, as well as the renormalization of the 
parameters g St y and F^ L ^ R . In the MS scheme the renormalization constants for the various 
fields can be extracted from the previous sections in this appendix. Splitting them again into 
class a) and class b) contributions, the non-vanishing ones read for the quarks 



5Z a = 5Z% L 



5Z b = 5Z b ? L 



5Z^ R = -^C F -, 
47r e 



5Z 



b,R 



•— C F - 

47r e 



and for the squarks 

5Z a = bZt 



kk 



<*s r 2 
— C F - 
47r e 



5Z b 



5Z kk> 



ryb a * n Z 

kk 4vr e 



a* 



K Xdy ~ X dk 



+ -X%x dl X$ 



-, (k + k') 



m d rr kd r k'd* , r kd r k'd* 
1 no tl ncl i 



DL L DR 



DR DL 



and finally for the gluino (Z 3 = 1 + 5Z§ + 8Z\) 



5Z^ = -^C A -, SZ b 

47T e 47r e 



a, 1 

ztrrif 



(119) 



(120) 



(121) 



The counterterm contribution to the vertex due to the renormalization of the fields can be 
written as 



sv: 



ct, fields 



ki 



= X - [bZ a + 5Z b + 5Z a + 5Z h + bZl + SZ^j V k 



tree 
ki 



(122) 



We now turn to the renormalization of the parameters g S) y and T DL / R matrices. For the 
former we write the connection between the bare- and renormalized version in the form 
g ba Y e = (1 + 5Zg a Y + $Z b s Y ) 9s,y while for the latter the corresponding connection reads 



r bare,ki _ \ ( M' , xrkk'\ r k'i 

1 DL/R - Z^\° ' D J * DL/R ■ 



(123) 



where 5T D k is an antihermitian 6x6 matrix. This procedure corresponds to a rotation of 
the squarks only. Therefore the renormalized T DL / R matrices still describe the transition 
between the super-CKM basis and the mass eigenstate basis. The counterterm contribution 
to the vertex due to the renormalization of the parameters then reads 



ryct,param 
ki 



(124) 
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irred j_ ^fields + ^param ^ bg ^ ^ ^ ^6 v and ^fefe' in t]|( , 



Requiring <5V™ + SV ki 
MS scheme to be (Z 9sX = 1 + SZ^ Y + SZ h gs Y 



SZl >Y = -^---(C A + C F ), 

Osl / 3 
4vr e V2 



SZ Ly = ±-AnC F + trn f ), 



(125) 



d *-D = 2 6Z kk' > 



where 5Z b k y H is given in (|120p . 
F.2 Renormalization of 



Stirred 



tree 
kk 



(126) 



We use the squark-squark-gluon-vertex to compute the one-loop renormalization constant for 
the strong coupling constant g St c of gauge type in the MS scheme. 

At the one- loop order, there are six one-particle irreducible diagrams contributing to the 
vertex correction. Two of them contain a gluino and a quark in the loop, one consists of 
a squark loop, while the remaining three feature a virtual gluon. For the vertex correction 
SV kk due to these irreducible diagrams we obtain (keeping only the UV-singularities) 

^c F -c A )\vSr + ^ic F \ 

where the first term on the r.h.s. is due to gluon corrections, while the second one results from 
gluinos and squarks. Additionally, there is a counterterm contribution &V kk which contains 
the following sources: renormalization of the squark and the gluino field, as well as the 
renormalization of g s ,G- Again, the renormalization constants for the various fields in the MS 
scheme can be extracted from the previous sections in this appendix. Unlike for the squark- 
quark-gluino-vertex, only the diagonal renormalization constants for the squark fields are 
involved in the squark-squark-gluon-vertex. The gluon-field renormalization constant reads 
(Z 3 = 1 + 5Z$ + bZ\) 



bZ% 



Air 



tr n 



1 



5Z\ 



4tt 



2 2 
-trn f + -C A 



1 



(127) 



where the first expression stems from gluon, quark and ghost contributions in the corre- 
sponding gluon self-energy, while the second one is due to gluino and squark loops. Note 
that the second term corresponds to the UV-singular part of bZ\ in (|116p . The diagonal 
renormalization constants for the squark fields are contained in ()120p . 

The counterterm contribution to the vertex due to the renormalization of the fields and 
of g S: c can be written as (Z gs 



svft 



,G 



6Z a + 5Z< 



l + 6Z« siG + 5Zl G ) 



+ l -5Zt+ l -5Z 



3 + $ Z g s ,G 



+ sz; 



g s ,G 



i/tree 
v kk ■ 



(128) 



Requiring SV kk ed + SV k l to be finite, fixes SZg sG and SZg gG in the MS scheme to be 



x 7 a 



2 11 ' 

- 3 tm f --C A 



6Z i,G 



47T 



1 C A 



(129) 
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